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Abstract 

We present a preliminary measurement of the time-dependent CP-violating asymmetry in decays 
of neutral B mesons to the final states D*^t:^, using approximately 178 million BB events recorded 
by the BABMl experiment at the PEP-II e^e~ storage ring. Events containing these decays are 
selected with a partial reconstruction technique, in which only the high-momentum vr^ from the B 
decay and the low-momentum vr^ from the D*^ decay are used. We measure the amplitude of the 
asymmetry to be —0.041 it 0.016 (stat.) it 0.010 (syst.) and determine bounds on | sin(2/3 + 7)|. 
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1 INTRODUCTION 



The Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing matrix [1] gives an explanation of CP 
violation and is under experimental investigation aimed at constraining its parameters. A crucial 
part of this program is the measurement of the angle 7 = arg (—yudYub/^cd^d)) unitarity 
triangle related to the CKM matrix. The decay modes D*^7r='= have been proposed for 

use in measurements of sin(2/3 + 7) [2], where /3 = arg{—V^V*f^/Vf.^V^D is well measured [3]. In 
the Standard Model the decays D*^tt~ and B^ D*^-k~ proceed through the b ucd 

and b cud amplitudes Au and Ac. The relative weak phase between the two amplitudes in the 
usual Wolfenstein convention [4] is 7. When combined with B^B^ mixing, this yields a weak phase 
difference of 2/3 + 7 between the interfering amplitudes. 
The decay rate distribution for B — > D*^ir^ is 

e-l^*l/^ r , 1 
P^{At) = — X Its'- sin(AmAt) t vC cos(AmAi) | , (1) 



where r is the B^ lifetime averaged over the two mass eigenstates, Am is the B^ — B^ mixing 
frequency, and At is the difference between the time of the B — *■ D*^Tr^ (B^ec) decay and the 
decay of the other B (i?tag) iii the event. The upper (lower) sign in Eq. (1) indicates the flavor of 
the Stag as a 5° (B^), while r? = +1 (-1) and C = + (-) for the Brec final state L»*"7r+ {D*+7r-). 
The parameters C and are given by 

C^'t^, 5±^-^^sin(2/3 + 7±fe*.). (2) 

Here 6d*tt is the strong phase difference between A^ and Ac and vd'-k = \Au/Ac\. Since A^ is 
doubly CKM-suppressed with respect to Ac, one expects rD*-jT ~ 0.02. 

We report a study of the CP-violating asymmetry in B^ D*^Tr^ decays using the technique 
of partial reconstruction, which allows us to analyze a large sample of signal events. We use 
approximately twice the integrated luminosity used in our previous analysis of this process [5] , and 
employ an improved method to eliminate a measurement bias, as described in Sec. 3.4.2. 



2 THE BABAR DETECTOR AND DATASET 

The data used in this analysis were recorded with the BABAR detector at the PEP-II storage ring, 
and consist of 165.6 fb~^ collected on the T'(4S) resonance (on-resonance sample), and 16 fb~^ 
collected at an e~^e~ CM energy approximately 40 MeV below the resonance peak (off-resonance 
sample). Samples of simulated Monte Carlo (MC) events with an equivalent luminosity 4 times 
larger than the data are analyzed through the same analysis procedure. 

The BABAR detector is described in detail in Ref. [8]. We provide a brief description of the 
main components and their use in this analysis. Charged-particle trajectories are measured by a 
combination of a five-layer silicon vertex tracker (SVT) and a 40-layer drift chamber (DCH) in a 
1.5 T solenoidal magnetic field. Tracks with low transverse momentum can be reconstructed in the 
SVT alone, thus extending the charged-particle detection down to transverse momenta of about 
50 MeV/c. We use a ring-imaging Cherenkov detector (DIRC) for charged-particle identification 
and augment it with energy-loss measurements from the SVT and DCH. Photons and electrons are 
detected in a CsI(Tl) electromagnetic calorimeter (EMC), with photon-energy resolution ge/E = 
0.023(£'/GeV)~^/^ 0.019. The instrumented flux return (IFR) is equipped with resistive plate 
chambers to identify muons. 
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3 ANALYSIS METHOD 



3.1 PARTIAL RECONSTRUCTION OF B« Z)*T7r± 

In the partial reconstruction of a B° — >^ D*^-!:"^ candidate (-Brec), only the hard (high- momentum) 
pion track iTh from the B decay and the soft (low-momentum) pion track vr^ from the decay 
D*~ — > -D°7r~ are used. Applying kinematic constraints consistent with the signal decay mode, we 
calculate the four-momentum of the unreconstructed, "missing" D, obtaining its flight direction to 
within a few degrees and its invariant mass mmiss [6]- Signal events peak in the mmiss distribution 
at the nominal mass Mjjo with an r.m.s. of 3 MeV/c^ (Fig- !)• 

3.2 BACKGROUNDS 

In addition to D*^tt^ events, the selected event sample contains the following kinds of 

events: B — > D*^ p^; BB background peaking in mmiss) composed of pairs of tracks coming from 
the same B meson, with the iTg originating from a charged D* decay, excluding B D*^p^ decays; 
combinatoric BB background, defined as all remaining BB background events; and continuum 
e~^e~ qq, where q represents a u, d, s, or c quark. We suppress the combinatoric background 
with selection criteria based on the event shape and the D* helicity-angle. We reject tt/j candidates 
that are identified as leptons or kaons. All candidates must satisfy the requirement 1.81 < mmiss < 
1.88 GeV/c^. Multiple candidates are found in 5% of the events. In these instances, only the 
candidate with the mmiss value closest to M^o is used. 

3.3 DECAY TIME MEASUREMENT AND FLAVOR TAGGING 

To perform this analysis, At and the flavor of the -Stag must be determined. We tag the flavor of 
the Stag using lepton or kaon candidates. The lepton CM momentum is required to be greater than 
1.1 GeV/c to suppress "cascade" leptons that originate from charm decays. If several flavor-tagging 
tracks are present in either the lepton or kaon tagging category, the only track of that category used 
for tagging is the one with the largest value of 6t, the CM angle between the track momentum 
and the missing D momentum. The tagging track must satisfy cos^t < Ct, where Ct = 0.75 
{Ct = 0.50) for leptons (kaons), to minimize the impact of tracks originating from the D decay. If 
both a lepton and a kaon satisfy this requirement, the event is tagged with the lepton only. 

We measure At using At = (zrec — -2tag)/(7/3c), where Zrec (-^tag) is the decay position of the 
Bj-ec (-Stag) along the beam axis (z) in the laboratory frame, and the e'^e~ boost parameter 7/? is 
continuously calculated from the beam energies. To find Zrec, we use the vr/i track parameters and 
errors, and the beam-spot position and size in the plane perpendicular to the beams (the x — y 
plane). We find the position of the point in space for which the sum of the contributions from 
the 7r?i track and the beam spot is a minimum. The z coordinate of this point determines z^^c- In 
lepton-tagged events, the same procedure, with the tt^ track replaced by the tagging lepton, is used 
to determine Ztag- 

In kaon-tagged events, we obtain ztag from a beam-spot-constrained vertex fit of all tracks in 
the event, excluding tt/j, tt^ and all tracks within 1 rad of the D momentum in the CM frame. If 
the contribution of any track to the of the vertex is more than 6, the track is removed and the 
fit is repeated until no track fails the < 6 requirement. 

The At error a/^t is calculated from the results of the z^ec and ztag vertex fits. 
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3.4 PROBABILITY DENSITY FUNCTION 



The analysis is carried out with a series of unbinned maximum-hkeUhood fits, performed simulta- 
neously on the on- and off-resonance data samples and independently for the lepton-tagged and 
kaon-tagged events. The probability density function (PDF) depends on the variables mmiss; ^t, 
o'At, F, Sj, and s^, where F is a Fisher discriminant formed from fifteen event-shape variables that 
provide discrimination against continuum events [6], = 1 (—1) when the -Btag is identified as a 

{B^), and = 1 (—1) for "unmixed" ("mixed") events. An event is labeled unmixed if the tt/j 
is a 7r~(7r+) and the Stag is a B^{B^), and mixed otherwise. 

The PDF for on-resonance data is a sum over the PDFs of the different event types: 

i 

where the index i = {D*-ir, D* p, peak, comb, qq} indicates one of the event types described above, 
fi is the relative fraction of events of type i in the data sample, and Pi is the PDF for these events. 
The PDF for off-resonance data is Pqq. The parameter values for Pj are different for each event 
type, unless indicated otherwise. Each Pi is a product. 

Pi = Miim^ss) HP) Tii^t, (7At, s„sj, (4) 

where the factors in Eq. (4) are described below. 

3.4.1 mmiss AND F PDFs 

The mmiss PDF for each event type i is the sum of a bifurcated Gaussian plus an ARGUS function: 

-Mi("T-miss) = fi Gi{mmiss) + (1 " /f ) A("imiss), (5) 

where is the bifurcated Gaussian fraction. The functions Qi and Ai are 
^^^"^^ \exp[-(x-M,)V2a|J , x > M, ' 



A{x) oc 



xJl - (x/M/^)^ exp 



(x/M^)' 



x<M/^ 



, x> 



(7) 



I 1 



where Mi is the peak of the bifurcated Gaussian, (jli and cjRi are its left and right widths, is 
the ARGUS exponent, Mf^ is its end point, and the proportionality constants are such that each 
of these functions is normalized to unit area. The mmiss PDF of each event type has different 
parameter values. 

The Fisher discriminant PDF Ti for each event type is parameterized as a bifurcated Gaussian, 
as in Eq. (6). The parameter values of Tb*-k-, ^D*p-, -^peak; and .Fcomb are identical. 

3.4.2 SIGNAL At PDFs 

The At PDF r^«^(At, cta*, s^, s^) for signal events corresponds to Eq. (1) with 0(r|)*^) terms 
neglected, and with additional parameters that account for several experimental effects, described 
below. 
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The first effect lias to do with the origin of the tagging track. In some of the events, the tagging 
track originates from the decay of the missing D. These events are labeled "missing-^) tags" and 
do not provide any information regarding the flavor of the -Stag- In lepton tag events we further 
distinguish between "direct" tags, in which the tagging lepton originates directly from the decay of 
the i?tag) and "cascade" tags, where the tagging lepton is a daughter of a charmed particle produced 
in the Stag decay. Cascade and direct tags have different mistag probabilities due to the different 
physical origin of the tagging track. In addition, the measured value of ztag in cascade-lepton tags is 
systematically larger than the true value, due to the finite lifetime of the charmed particle and the 
boosted CM frame. This creates a correlation between the tag and vertex measurements that we 
address by considering cascade-lepton tags separately in the PDF ^. In kaon tags, ztag is determined 
using all available 5tag tracks. Therefore, the effect of the tagging track on the ztag measurement 
is small, and no distinction between cascade and direct kaon tags is needed. 

The second experimental effect is the finite detector resolution in the measurement of At. We 
address this by convoluting the distribution of the true decay time difference Attr with a detector 
resolution function. Putting these two effects together, the At PDF of signal events is 

r^.^(At,aAt,s„s^) = {l + AeD*nS,)^fj,*^J dAtt,T'j,,^{Attr,s„sjn^,^{At - Atu,aM), 

j 

(8) 

where 2AeD*-K is the relative difference between the detection efficiencies of positive and negative 
leptons or kaons, the index j = {dir, cas, miss} indicates direct, cascade, and missing-D tags, and 
/|)«^ is the fraction of signal events of tag-type j in the sample. We set f^l^ = 1 — f^S^ — /™*^- For 
kaon tags fjf*^ = 0. The function T£,*^(Attr; "^m) is the Attr distribution of tag-type j events, 
and 7^£),^(At — Attrj ^"Af ) is their resolution function, which parameterizes both the finite detector 
resolution and systematic offsets in the measurement of Az, such as those due to the origin of the 
tagging particle. The parameterization of the resolution function is described in Sec. 3.4.4. 
The Attr PDF for missing-D tags is 

rg'i-lAttr, sj = ' ^ {1 + (1 - 2pn^,)} , (9) 

where Pd^tt is the probability that the charge of the tagging track is such that it results in a mixed 
flavor measurement. 

The functional form of the direct and cascade tag Atn PDFs is 

l^*tr|/T£,* 

_J 1 _ o A, , 



Tj,,^{Attr,s„sJ = — I l-s,A, 

(1 - 2a;|j,^) cos(AmD*7rAttr) 
-S^D*n sin(Am£)*7rAttr)|, (10) 

where j = {dir, cas}, the mistag rate ui'jj*^ is the probability to misidentify the flavor of the Stag 
averaged over and and Aw^*^ is the B'^ mistag rate minus the B^ mistag rate. The factor 
Sj^t^ describes CP violation due to interference between b — > ucd and b — cud amplitudes in both 
the Bj-ec and the Stag decays: 



^In Ref. [5] we corrected for the bias caused by this effect and included a systematic error due to its uncertainty. 
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where a^)*^, bD*-;^, and cd*-,^ are related to the physical parameters through 

aD*7v = 2rD*nSin{2P + -f) cos 5d*tv, 
bo^Tv = 2r'sin(2/3 + 7)cosy, 

cd*tv = 2cos(2/3 + 7)(rD*7rSin(5D*7r-''^'sin(5'), (12) 

and r' {5') is the effective magnitude of the ratio of amplitudes (strong phase difference) between 
the b ucd and b cud amplitudes in the i^tag decay. This parameterization is good to first 
order in r£)*7r and v'. Xhe inclusion of v' and S' in the formalism accounts for cases where the 
-Stag undergoes a b ^ ucd decay, and the kaon produced in the subsequent charm decay is used for 
tagging [7]. In lepton-tagged events r' = (and hence bD*Tr = 0). 



3.4.3 BACKGROUND At PDFs 

The At PDF of B ^ D*^ has the same functional form and parameter values as the signal PDF, 
except that the CP parameters aD*pi bD*p, and CD*p are set to and are later varied to evaluate 
systematic uncertainties. The validity of the use of the same non-CP parameters for T'j^tp and T'j-,t.^ 
is established using simulated events, and stems from the fact that the tt/j momentum spectrum in 
the B D*^p^ events that pass our selection criteria is almost identical to the signal spectrum. 
The At PDF of the peaking background accounts separately for charged and neutral B decays: 

Tpeak{^t,aAt,S„sJ = (l + S^Aepeak) 

X {r°'peak + / dAtu rp+,k(Mr, S„ sj 7^+ 3^(Ai - Aitr, CTAt)] , (13) 

where T'^pcak the functional form of Eq. (8) and the subsequent expressions, Eqs. (9-12), 
but with all D*7r-subscripted parameters replaced with their peak-subscripted counterparts. The 
integral in Eq. (13) accounts for the contribution of charged B decays to the peaking background, 
with 

-|Attr|/r+ . 

e peak / \ 

r+,k(Attr, sj = -—^ (l - Au;+ , (14) 

peak 

and Tl^^^y^{At — Att^, a At) being the 3-Gaussian resolution function for these events. 

Convergence of the fit with the parameterization of T^*^ relies on our ability to independently 
determine fj^l^^ using the angle between the tTs and the D*, as described later in this note. That 
determination may not be carried out reliably for the combinatoric BB background, due uncer- 
tainties in simulating the random nature of the reconstructed track combinations. As a result, the 
four parameters /comb' '^comb' ^'^comb' ^^'^ Pcomb are not independent, and can be reduced to the 
set of three parameters 



/ _ dir /-i _ j-dir \ , J comb 

'^comb "^comb \ J comb J ' 

dir 



A'^comb - A^comb (1 " /comb)' 

■^comb = /comb(-'- ~ 2 Pcomb)- (1^) 

With these parameters and f^omb ~ ^' combinatoric BB background At PDF becomes 

^comb(Mo-At,Sj,v) = (1 + St A^comb) j dAtu Tcombi^kv, S^, sj Tlcombi^t - ^hr, CTAt) , (16) 
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where 7^comb(^* ~ ^hr^f^At) is the 3-Gaussian resolution fucntion and 

g— |Attr|/Tcomb f 

rcomb(Attr, St> ■Sm) = 7 \ 1 " ^^^comb + Sm^comb 

^"^comb I 

+Sm (1 - 2Wcomb) COs(AmcombAitr) 



-'^comb sin(AmcombAitr)|, (17) 



with 

'5comb = (1 - 2a;comb) (StOcomb + SmCcomb) + S^S^bcomhi^ " ^t^'^comb)- (18) 

As in the case of Td*p, the CP parameters apeak, ^'peak, and Cpeak are set to and are later varied 
to evaluate systematic uncertainties. Parameters labeled with superscripts "peak" or "comb" are 
empirical and thus do not necessarily correspond to physical parameters and may have different 
values from those of the D*7r-labeled parameters. 

The PDF Tqq for the continuum background is the sum of two components, one with a finite 
lifetime and one with zero lifetime: 

T^^(At, aAt, sj = (1 + Ae,^) J dAttr T,^(Attr, s,, s^) 7^,^(A^ - Att„ (Ta*), (19) 



with 

TqqiAtu, = (1 - 4)^- (1 - Augg) + f^qS{Atu), (20) 

where /^^ is the fraction of zero-lifetime events. 

3.4.4 RESOLUTION FUNCTION PARAMETERIZATION 

The resolution function for events of type i and optional secondary- type j {j = {dir, cas, miss} 
for lepton-tagged signal events and j = {+, 0} for the BB background types) is parameterized as 
the sum of three Gaussians: 

niiAt - Attr, aAt) = n Q'^liU, CTAt) + (1 - n - rl) Q'^Htr^aAt) + n Oi^tr, aAt), (21) 

where tr is the residual At — Attr, and G'^l, and Q°l are the "narrow", "wide", and "outlier" 
Gaussians. The narrow and wide Gaussians have the form 



Q^l{tr,aAt) = exp 

V27r s^l a At 



^ _^r-y^laA^ 

^ Hsf^laAt)' ^ 



(22) 



where the index k takes the values k = n,w for the narrow and wide Gaussians, and 6*^j and s'^j 
are parameters determined by fits, as described in Sec. 3.5. The outlier Gaussian has the form 



1 ( i^r-b^i) 

G°i {tr, aAt) = ^ exp - ^ ^ 



27rs 



V 



2(.°])2 



2^ 

(23) 



where in all nominal fits the values of h°l and s°l are fixed to ps and 8 ps, respectively, and are 
later varied to evaluate systematic errors. 



13 



3.5 ANALYSIS PROCEDURE 



The analysis proceeds in four steps involving unbinned maximum likelihood fits to the data: 

1. In the first step, we determine the parameters fo'p + fo'iv^ /peak) a-i^d /comb of Eq. (3). In 
order to reduce the reliance on the simulation, we also obtain in the same fit the parameters 

of Eq. (5), eqq of Eq. (7), ctl for the signal (Eq. 6) and all the parameters of Tqq, and J^^q 
(the latter applies to all BB event types). This is done by fitting the data with the PDF 

Pi = Mi(rnr^,s)J'i{F), (24) 

instead of Eq. (4), i.e. by ignoring the time dependence. The fraction fqq of continuum events 
is determined from the event yield of the off-resonance sample and its integrated luminosity 
relative to the on-resonance sample. All other parameters of the Mi PDFs and the value of 
fD*-K/{fD*'K + fD*p) are obtained from the MC simulation. 

2. In the second step, we repeat the fit of the first step for data events with cos > Cy, to obtain 
the fraction of signal events in that sample. Given this fraction and the relative efficiencies 
for direct, cascade, and missing- £) signal events to satisfy the cos^t < Ct requirement, we 
calculate /j^i^- We also calculate the value of po*n from the fractions of mixed and unmixed 
signal events in the cos 9t > Ct sample relative to the cos 9t < Ct sample. 

3. In the third step, we fit the datci events in the sidebcinci 1.81 <C ^miss 

< 1.84 GeV/c2 with the 

3-dimensional PDFs of Eq. (4). The parameters of Mi{mraiss) and J^i{F), and the fractions 
fi are fixed to the values obtained in the first step. From this fit we obtain the parameters 
of Tp'omb) as well as those of T^^. 

4. In the fourth step, we fix all the parameter values obtained in the previous steps and fit 
the events in the signal region 1.845 < "mmiss < 1-880 GeV/c^, determining the parameters 
of T[)*^ and T^^. Simulation studies show that the parameters of T^^^^y^ are independent of 
fJ^miss) enabling us to obtain them in the sidebad fit (step 3) and then use them in the signal- 
region fit. The same is not true of the T^^ parameters; hence they are free parameters in the 
signal-region fit of the last step. 



4 RESULTS 

The fit of step 1 finds 16060 it 210 signal 5^ D*^Tr^ events in the lepton-tag category and 
57480 ± 540 in the kaon-tag category. The mmiss and F distributions for data are shown in Figs. 1 
and 2, with the PDFs overlayed. 

The results of the signal region fit for lepton-tagged events are summarized in Table 1, and 
the plots of the At distributions for the data are shown in Fig. 3. Results of the fit for the kaon- 
tagged events are shown in Table 2 and Fig. 4. For each of the plots in this figure we calculate the 
Kolmogorov-Smirnov probabilities for the PDF and data to originate from the same distribution, 
as a way to verify the goodness of the fit. The probabilities are 41%, 99%, 99%, and 26%. Fig. 5 
shows the raw, time-dependent CP asymmetry 

Ns -i(At)-Ns —i(At) 

In the absence of background and with high statistics, perfect tagging, and perfect At measurement, 
A{At) would be a sinusoidal oscillation with amplitude aD*TT- 
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Figure 1: The m-miss distributions for on-resonance lepton-tagged (left) and kaon-tagged (right) 
data. The curves show, from bottom to top, the cumulative contributions of the continuum, 
peaking BB, combinatoric BB, B D*^p^, and B^ D*^Tr^ PDF components. 




Figure 2: The F distributions for on-resonance lepton-tagged (left) and kaon-tagged (right) data. 
The contributions of the BB (dashed line) and the continuum (dot line) PDF components are 
overlayed, peaking at approximately —0.6 and —0.1, respectively. Also overlayed is the total PDF. 
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Parameter 


Value 


signal 


2rD*7r sin(2/3 + 7) cosSd^w 
2rD*n cos(2/3 + 7) sin SD*n 




-0.048 ± 0.022 
-0.015 ±0.036 


— mixing frequency 

B^ lifetime 
mistag (direct tags) 


TD*TT 


0.537 ±0.011 ps-i 
1.435 ±0.019 ps 
0.008 ± 0.0.003 


bias of g^i^ 

bias of 

bias oigi, 
fraction of ^Jir 
fraction of ^^^^ 
width of g^ir 

width of gj,^,. 


un cas 
un dir 

'^D*7T 

uw dir 

fn dir 
JD*Tr 
to dir 

J D*TT 

„n dir 

^D*Tr 
„w dir 


-0.43 ± 0.22 
0.000 ± 0.034 

0.04 ± 0.40 
0.918 ± 0.065 
0.005 ±0.005 
1.047 ± 0.063 

2.43 ± 0.68 


coulinuuni 


qq effective lifetime 
mistag 
mistag difference 
fraction of 0-lifetime events 


Tqg 

AUJqg 

Jqq 


1.12 ± 0.22 ps 
0.343 ± 0.011 
-0.095± 0.022 
0.775 ± 0.037 


bias of 
bias of g'"" 
fraction of g'^ 
fraction of g° 

width (>[ g" 


6"- 

qq 
qq 

j!n 
Jqq 

f - 

qq 

•s"- 


0.047 ± 0.054 
-1.73± 0.83 
0.848 ± 0.055 
0.072 ± 0.015 

0.982 ± 0.052 



Table 1: Results of the fit of the lepton-tagged events in the signal region 1.845 < mmiss < 
1.880 GeV/c^. Errors are statistical only. 
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Figure 3: At distributions for the lepton-tagged events separated according to the tagged flavor of 
-Stag and whether they were found to be mixed or unmixed: a) unmixed b) mixed c) B^ 
unmixed, d) B^ mixed. The curves show the PDF, calculated with the parameters obtained by 
the fit. Also shown is the PDF of the total background. The requirements mmiss > 1-855 GeV/c^, 
F < are applied in order to reduce the background. 
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Parameter 


Value 


signal 


2rD'-TT sin(2/3 + 7) cos 5d*tt 




-0.033 ± 0.023 


2r' sin(2/3 + 7) cos 5' 




-0.004 ±0.012 


2cos(2/3 + 7)(r£)*7r smSD*-iT — r' sind') 




0.019 ± 0.023 


— mixing frequency 




0.4716 ± 0.0087 ps"^ 


B^ effective lifetime 


TD*n 


1.379 ± 0.012 ps 


mistag 




0.226 ± 0.0036 


mistag difference 




-0.0301 ± 0.0052 


efficiency difference 


AeD*7r 


-0.0126 ± 0.0055 


bias of 


"D*tt 


-0.287 ±0.022 


bias of g"" 


"D*tt 


0.044 ± 0.067 


fraction of ^" 


fn 

JD*Tr 


0.909 ± 0.040 


fraction of g° 


/D*7r 


0.014 ± 0.002 


width of 


„n 

^D*w 


1.109 ±0.039 


width of g"" 


*D*7r 


0.01 ± 0.23 


continuum 


qq effective lifetime 


Tqq 


0.546 ± 0.037 ps 


mistag of non-O-lifetime events 




0.000 ± 0.002 


mistag of 0-lifetime events 


4 


0.334 ±0.010 


fraction of 0-lifetime events 


Jqq 


0.768 ± 0.021 


bias of 


6"- 

qq 


0.013 ± 0.007 


bias of g"" 


qq 

Jqq 

f - 

Jqq 


0.130± 0.047 


fraction of g^ 


0.853 ± 0.026 


fraction of g° 


0.013 ± 0.001 


width of 


s"- 

</</ 


1.029 ± 0.013 


width of g"" 


qq 


1.90 ± 0.11 



Table 2: Results of the fit of the kaon-tagged events in the signal region 1.845 < mmiss < 
1.880 GeV/c^. Errors are statistical only. 



18 



in 


. 1 1 1 . 1 . , . 1 1 1 , 1 1 1 1 


,1.1.,.. 


Tr3500 

CM 






O 




^3000 
c 


-BaBar |l 




a> 






l"2500 


-preliminary jM 




2000 






1500 






1000 






500 








-10 -5 5 10 

A t (ps) 




A t (ps) 




Figure 4: At distributions for the kaon-tagged events separated according to the tagged flavor of 
Stag and whether they were found to be mixed or unmixed: a) unmixed b) mixed c) B'^ 
unmixed, d) B^ mixed. The curves show the PDF, calculated with the parameters obtained by 
the fit. Also shown is the PDF of the total background. The requirements mmiss > 1.855 GeV/c^, 
F < are applied in order to reduce the background. 
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Figure 5: Raw asymmetry for (a) lepton-tagged and (b) kaon-tagged events. The curve represents 
the projection of the PDF for the raw asymmetry. The requirements m^iss > 1.855 GeV/c^, F <0 
are apphed in order to reduce the background. A non-zero value of aD*n would show up as a 
sinusoidal asymmetry, up to resolution and background effects. The offset from the horizontal axis 
in the kaon-tag plot is due to the nonzero value of A€D*n- 
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bource 


Error (xlO ^) 








1. Step 1 fit 


0.04 


0.04 


2. Sideband statistics 


0.08 


0.08 




0.02 


0.02 


4. p 


0.02 


0.02 


5. MC statistics 


0.6 


1.2 


6. Beam spot size 


0.10 


0.10 


7. Detector z scale 


0.03 


0.03 


8. Detector alignment 


0.4 


0.8 


9. Combmatoric background CP content 


0.25 


0.22 


10. Peaking background CP content 


0.36 


0.38 


11. D*p CP content 


0.53 


0.52 


12. Peaking background 


0.21 


0.31 


lo. Signal region/ side band diiierence 


0.0003 


0.002 


14. B{B D*^p^) 


0.17 


0.33 


15. Variation of r/j*^ and Ami3*7r 


0.21 


0.95 


Total systematic error 


1.04 


1.89 


Statistical uncertainty 


2.2 


3.6 



Table 3: Systematic errors in af),^ and cfj^*^ for lepton-tagged events. 

5 SYSTEMATIC STUDIES 

The systematic errors are summarized in Tables 3 and 4 for lepton- and kaon-tagged events, re- 
spectively. Each item below corresponds to the item with the same number in Tables 3 and 4. 

1. The statistical errors from the fit in Step 1 are propagated to the final fit, taking their 
correlations into account. It also includes the systematic errors due to possible differences 
between the PDF line shape and the data points in the kinematical fit. 

2. The statistical errors from the mmiss sideband fit (Step 3) are propagated to the final fit 
(Step 4), taking their correlations into account. 

3-4. The statistical errors from the Step 2 fits are propagated to the final fit. 

5. The statistical errors associated with all the parameters obtained from MC are propagated 
to the final fit. 

6. The effect of uncertainties in the beam-spot size on the vertex constraint is estimated by 
increasing the beam spot size by 50 //m. 

7. The effect of the uncertainty in the measured length of the detector in the z direction is 
evaluated by applying a 0.6% variation to the measured values of At and (T/\f 

8. To evaluate the effect of possible misalignments in the SVT, signal MC events are recon- 
structed with different alignment parameters, and the analysis is repeated. 
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Source 


Error (xlO 










K 


1. Step 1 fit 


0.10 


0.04 


0.04 


2. Sideband statistics 


0.40 


0.12 


0.44 


Q J? miss 
JD*tv 


0.02 


negl. 


negl. 


4. p 


0.02 


negl. 


negl. 


5. MC statistics 


0.8 


0.4 


0.9 


6. Beam spot size 


0.07 


0.13 


0.06 


7. Detector z scale 


0.02 


negl. 


0.03 


8. Detector alignment 


0.41 


0.14 


0.74 


9. Combinatoric background CP content 


0.80 


0.56 


0.72 


10. Peaking background CP content 


0.29 


0.17 


0.27 


11. D*p CP content 


0.57 


0.58 


0.58 


12. Peaking background 


0.21 


0.41 


0.31 


lo. Signal region/ sideband dinerence 


0.04 


0.03 


0.05 


14. ^(5 ^ 


0.17 


0.22 


0.33 


15. Variation of td*tt and AmD*Tr 


0.26 


0.16 


0.05 


Total systematic error 


1.47 


1.06 


1.63 


Statistical uncertainty 


2.3 


1.2 


2.3 



Table 4: Systematic errors in a^,^, ^D'n^ and c^.^ for kaon-tagged events. 



9-11. The CP parameters of the B D*^p^, peaking, and combinatoric BB background are fixed 
to in the fits. To study the effect of possible CP violation in these backgrounds, their CP 
parameters are varied in the range ±0.04 and the step-4 fit is repeated. 

12. The uncertainty due to the parameters of T^^^^^ is evaluated by fitting the simulated sample, 
setting the parameters of T^^^y. to be identical to those of T^.^^^. 

13. The uncertainty due to possible differences between the At distributions for the combinatoric 
background in the mmiss sideband and signal region is evaluated by comparing the results of 
fitting the simulated sample with the T^^^ parameters taken from the sideband or the signal 
region. 

14. The ratio fD*p/ fD*-K is varied by the uncertainty in the corresponding ratio of branching 
fractions, obtained from Ref. [10]. 

15. The lifetime and mixing parameters from the fits are not very consistent with the world 
average values when only the statistical uncertainties from the fit are considered. However, 
the inconsistency is significantly reduced when the estimated systematic uncertainties on 
the td*-k and AmD*T^ parameters are included. To determine the systematic errors on the 
CP parameters due to td*tt and AmD*Tr, we repeat the fit with these parameters fixed to 
their world-average values from Ref. [10]. Because of the small correlations between these 
parameters and the CP parameters (< 2% for the kaon-tagged sample and < 9% for the 
lepton-tagged sample), the resulting changes in the CP parameters are small. These changes 
are assigned as systematic uncertainties. 
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6 PHYSICS RESULTS 



Summarizing the values and uncertainties of the CP parameters, we obtain the following results 
from the lepton-tagged sample: 

ai,.^ = 2r£,*^sin(2/3 + 7)cos(5D*,r = -0.048 ± 0.022 ± 0.010, 

ci,.^ = 2r£i*^cos(2/3 + 7)sin(5D*7r = -0.015 ± 0.036 ± 0.019. (26) 

The results from the kaon-tagged sample fits are 

a^*^ = 2rD*nSm{2j3 + -i)cos5D*n = -0.033 ± 0.023 ± 0.015, 
6^,^ = 2r' sin(2/3 + 7) cos y = -0.004 ± 0.012 ± 0.011, 
c^.^ = 2cos(2/3 + 7)(r£,.7rSin(5£,*^ -r'sinJ') = +0.019 ± 0.023 ± 0.016. (27) 

Combining the results for lepton and kaon tags gives the amplitude of the time-dependent CP 
asymmetry, 

ao-TT = 2r£i*,r sin(2/3 + 7) cos Sd'^k = -0.041 ± 0.016 {stat.) ± 0.010 {syst.), (28) 

where the systematic error takes into account correlations between the individual results. This 
result deviates from zero by 2.2 standard deviations. 

We use two methods for interpreting our results in terms of constraints on | sin(2/3 + 7)|. Both 
methods involve minimizing a function that is symmetric under the exchange sin(2/3 + 7) — >■ 
— sin(2/3 + 7), and applying the method of Ref. [11]. In the first method we make no assumption 
regarding the value of td*t^. For different values of rD*i:-, we vary 5d*i: and sin(2/3 + 7) so as to 
minimize the function 

3 

x2(r-D*^,sin(2/3 + 7),(5D*^) = ^ Ax.F^^iAxfc, (29) 

where Axj refers to the difference between the result of our measurement of a^*^, a^£)*^, or c^j*^ 
(Eqs. (27) and (26)) and the theoretical expressions given by Eq. (12). The measurements of h^*^ 
and c^*^ are not used in the fit, since they depend on the unknown values of r' and 5' . The 
measurement error matrix V is nearly diagonal, and accounts for correlations between the mea- 
surements due to correlated statistical and systematic uncertainties. The parameters determined 
by this fit are sin(2/3 + 7), which is limited to lie in the range [—1, 1], and 5d*-k- We then generate 
many parameterized MC experiments with the same sensitivity as reported here for different values 
of sin(2/3 + 7) and with 5d*-k = 0, which yields the most conservative lower limit on | sin(2/'3 + 7)]. 
The fraction of these experiments in which x^(sin(2/3 + 7)) — Xmin smaller than in the data is 
interpreted as the confidence level (CL) of the limit on | sin(2/9 + 7)|. The resulting 90% CL lower 
limit on | sin(2/3 -|- 7)! is shown as a function of td'tt in Fig. 6. This limit is always based on the 
more conservative of the two possibilities implied by the ambiguity | sin(2/3 + 7)| | cos 5 d*i:\- 

The second method assumes that r/5*7r can be estimated from the Cabibbo angle, the ratio 
of branching fractions B{B^ D*f7T^)/B{B^ L'*^7r+) [12], and the ratio of decay constants 
Id*/ Id* [13], yielding the measured value 

rgS- = 0.015ingt (30) 
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Figure 6: Lower limit on | sin(2/3 + 7)] at 90% CL as a function of ro^-w- 



This value includes our recent measurement of the branching fraction B{Df — > (pn^) [14]. In 
addition to the above experimental errors, we attribute a non-Gaussian 30% relative error to 
the theoretical assumptions involved in obtaining this value. To obtain the limits with these 
assumptions, we minimize the function 



(31) 



where the term A^(r£)*7r) takes into account both the Gaussian experimental errors of Eq. (30) and 
the 30% theoretical uncertainty [15]: 



rD*n — '-■01 D 

mm 



Q J ^meas \ 2 



OHDT 



e.,,. > 0.3 

ir^^_ < -0.3 



(32) 



where Cr^,*^ = {td'tt — r^t^) / rjji'^ . The parameters 

sin(2/3 + 7) = 0.97 ±1.22, 
6D'n = 0.24 ±1.23, 
rD*n = 0.020 ± 0.004 



(33) 



are determined in this fit. Due to the fact that the minimum value occurs close to the boundary 
of the physical region (| sin(2/3 + 7)| = 1), the errors in Eq. 33 are not relevant and in order to give 
a probabilistic interpretation the method of Ref. [11] has been used. The resulting confidence level 
as a function of the lower limit on | sin(2/3±7)|, evaluated using parameterized MC experiments as 
in the first method, is shown in Fig. 7. In particular, we find the limits | sin(2/3 ± 7)| > 0.75 (0.58) 
at 68% (90%) CL. The implied contours of constant probability for the apex of the unitary triangle 
appear in Fig. 8. 
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Figure 7: Confidence level as a function of the lower limit on | sin(2/3 + 7)| given r^i*,^ from Eq. (30). 



7 SUMMARY 



We present a preliminary measurement of the time-dependent CP asymmetry and parameters 
related to sin(2/3 + 7) in a sample of partially reconstructed — > D*~^Tr~ events. In particular, 
the amplitude of the measured asymmetry is 

a^,*^ = 2ri5.,rSin(2/3 + 7) cos (5z)*^ = -0.041 ±0.016 (stat.) ± 0.010 (syst.). (34) 

We interpret our results in terms of the lower limits | sin(2/3 + 7)] > 0.75 (0.58) at 68% (90%) 
CL, and extract limits as a function of the ratio r^.^r between the b — > ucd and b cud decay 
amplitudes. 
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Figure 8: Contours of constant probability (color-coded in percent) for the position of the apex of 
the unitary triangle to be inside the contour, based on the results of Fig. 7. The one- and two- 
standard deviation ranges of the world average sin 2/3 measurement (sin 2/3 = 0.736 it 0.049) are 
shown as gray lines. The cross represents the value and errors of the apex of the unitarity triangle 
from the standard CKMFitter fit [16]. 
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